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ABSTRACT 

We report on evolutionary calculations of the onset of mass transfer in AM CVn binaries, treating 
the donor's evolution in detail. We show that during the early contact phase, while the mass transfer 
rate, M, is increasing, gravity wave (GW) emission continues to drive the binary to shorter orbital 
period, P or b- We argue that the phase where M > and v > (y = l/P or b) can last between 10 3 and 
10 6 yrs, significantly longer than previously estimated. These results are applied to RX J0806+1527 
(P orb = 321 s) and RX J914+2456 (P orb = 569 s), both of which have measured v > 0. Thus, a v > 
does not select between the unipolar inductor and accretion driven models proposed as the source of 
X-rays in these systems. For the accretion model, we predict for RX J0806 that v ps 1.0 — 1.5 x ICC 28 
Hz s -2 and argue that timing observations can probe v at this level with a total w 20 yr baseline. We 
also place constraints on each system's initial parameters given current observational data. 
Subject headings: binaries: close — gravitational waves — stars:individual (RX J0806+1527, RX 
J1914+2456)— white dwarfs 



1. INTRODUCTION 

The true nature of the two candidate ultracom- 
pact binary system s originally discove r ed wi th ROSAT, 
RX J0806+ 1527 llBeuermann et alJ Il999() and RX 
J1914+2456 (|Haberl fc Motchlll995|) . has been a source 
of much controversy. Both systems are soft X-ray sources 
whose light curves are modulated on periods of 321 
s (RX J0806, her eafter J0806) and 569 s (RX J1914 
hereafter J19141 llMotch et all I19961: llsrael et alJ 119991: 
iBurwitz fc Reinschl 1200 ljk In both cases, the X-ray 
light is 100% modu l ated for roughly half the period 
{Cropper et alJll998t llsrael et al]ll999ft and the optical 
light is modula t ed on the same period as the X-ray light 
l|Ramsav et al.| |2000t | Israel et all l2002t iRamsav et al.l 
I2002albl : llsrael et al.ll2003l) with little evid ence for other 
periodicities (but see IRamsav et al.ll2006j) . It is widely 
believed, based partially on the stability and singularity 
of these periods, that the modulations ar e on the sys- 
tems ' orbital periods, P or b (although see iNorton et alJ 
120041 for an alternate interpretation). This would make 
J0806 and J1914 the two shortest period binaries known. 

The X-ray production mechanism in these systems 
has been much debated. Two competing models have 
come to the fo re: the unipolar inductor (UI) model 
( Wtj^et^aTJ 2002) a nd the direct-impact accretion model 
( Marsh fc Steeghsl I2002D . Both model these systems 
as white dwarf-white dwarf (WD-WD) binaries with 
Porb equal to the X-ray period. In the UI model, a 
magnetized, more massive primary WD spinning asyn- 
chronously with respect to the orbit induces an electric 
field in the secondary WD and drives a current between 
the two. X-rays are produced by resistive dissipation 
in the primary's atmosphere. In the accretion model, 
the secondary fills it Roche Lobe (Rl) and mass trans- 
fer occurs. The binary's compact geometry leads to the 
accretion stream directly impacting the primary and no 
accretion disk forms. Both models lead to a spatially 
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small X-ray production site on the primary, explaining, 
in p rinciple, the observ ed X-ray modulation (see, how- 
ever iBarros et al.ll2005T> . 

One potential means of distinguishing between the two 
models is the secular time derivative of P or b, or equiv- 
alcntly of the orbital frequency, v = l/P or b, v. The 
UI model predicts v > 0. Assuming fully degenerate 
donors and enforcing exactly that the donor's radius, R2 
equals its Roche radius: Rl = R2, the accretion model 
woul d predict v < 0. Determinations of v from bot h op- 
tical l)Hakala et al.ll2003l 120041 IRamsav et alJl2005D and 
X-ray observations (|Strohmaverll2002t 120031 12004 12005(1 
all show that v > over the entire epoch of observa- 
tions. On the assumption that these results reflect the 
secular z>, this has been taken as evidence against ac- 
cretio n models fe.g.JIsrael et al.1120031 LStrohmaverll2002L 
120041 IRamsav et al.ll2005lL However, the R 2 = R L con- 
straint requires a very large mass transfer rate M k, 
10~ 7 — 10 _5 M© yr -1 and if M is significantly below this 
"equilibrium" value, observing a v > is possible. The 
M could be lower than expected either due to a non- 
secular mechanism (|Marsh fc Nelemansll2005l) or during 
the m ass transfer turn-on phase IjWillems fc Kaloeeral 
|2005|) . However, the timescales on which the binary 
(re)establishes its equilibr ium M are estimated to be 
very shor t: 100-1000 vrs llMarsh fc , Nelemansl l2005|l or 
2-20 yrs ijWillems fc Kalogerall2005|) . making observing 
such phases unlikely. 

Here we reexamine the suggestion J0806 and J1914 
are direct-impact accretors seen during the M turn-on 
phase in the context of stellar evolution models devel- 
oped to realistically address the complete evolution of 
AM CVn binaries. In §2, we show that the population 
of AM C Vns forming through a double-degenerate chan- 
nel fsee iNelemans et al.ll200ll iDelove et al.ll2005|) nat- 
urally produces systems whose contact phase evolution 
give Porb and v consistent with both J0806 and J1914. 
ccl Thc maximum M is sufficiently high in these systems 



2 



Dcloyc & Taam 



that Rl must penetrate below the donor's photosphere, 
naturally increasing the M growth times to s» 10 2 — 10 6 
yrs. We explore in §3 the range of initial conditions that 
will lead to J0806 and J1914. We determine constraints 
on these initial conditions and present predictions for v 
in J0806 as a diagnostic tool. Finally, we present our 
discussion in §4. 

2. MASS TRANSFER TURN-ON IN AM CVN SYSTEMS 

We take as our model for J0806 and J1914 a WD- 
WD binary with primary and secondary (donor) masses 
Mi and M2, where the donor is pure He. Binary evolu- 
tion through two common envelope (CE) phases can pro- 
duce such systems that gravity wave (G W) emission can- 
drive into contact within a Hubble time ijNelemans et alJ 
120011) . During this pre-contact phase, the donor cools 
and contracts. Variations in the donor's pre-CE entropy 
and time to contact lead to a range in the donor's de- 
generacy and radius, P 2 , at contact ijDelove fe Bildstenl 
120031 IDelove et alJl2005(l . For fixed M 2 , less degenerate 
donors have larger R2, leading to longer P or b at contact, 

Porb.i ■ 

We parameterize the system's state at contact by M^i, 
M 2 ,i, and i?2,i- The range of Mi j , M? j pa irs we con- 
sider are taken from iNelemans et al.l (|2001|) . The i?2,i 
range is dete r mined using the prescription outlined in 
IDelove et al.l 1)20051) with slight modifications detailed 
in a companion paper (Deloye et al. 2006, in prepara- 
tion). For this population, fully degenerate donors pro- 
duce P rb,i < 6 min; a realistic treatment of pre-contact 
donor cooling gives P or b,i extending well above 10 min 
ijDelove et al.ll200"> :. Thus an accretion model can easily 
accommodate the P or b of both J0806 and J1914 as early 
contact systems. 

We have carried out a range of evolutionary calcula- 
tions beginning from the initial conditions determined 
above. We utilized a coupled stellar/binary evolution 
code specifically developed for these calculations and fol- 
lowed each model's evolution from pre-contact through 
the late stages of the mass transfer phase. We assume 
conservative mass transfer, circular orbits, and that GW 
emission is the only source of orbita l angular mom entum, 
J, loss. We use the prescription of lRitterl l)1988|) to cal- 
culate M when Rr, > R2 and adopt the prescription of 
iKolb fc Rlttel IpMf when R L < R 2 . 

Under these assumptions, v is given by 



where j/J is given in lLandau fc Lifshitzl l)1971|) . q = 

M 2 /M 1 , and M 2 = -M < 0. As j/J < 0, the first term 
in equation JJJ drives v > 0; a sufficiently high M is 
required for i> < 0. After the donor comes into contact, 
there is a period of time in which both M2 7^ and v > 0. 
How quickly i> reaches depends on the growth rate of 
M. 

To place J0806 and J1914 in the context of the turn- 
on phase in AM CVn binaries, we show in Figure ^ the 
v-Poth evolution for a sample of our calculations. Solid 
lines originating on the same dotted line share the same 
Mi,i, M24, but different i?2,i- The dotted lines show 
the corresponding relation for GW losses alone. The di- 
amonds show the measured v for J0806 (3.63 x 10" 16 
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Fig. 1. — Comparison between the predicted, secular evolution of 
v vs. P nr h and the measu red values of i> in RX J0806 and RX J1914 
IStrohm avcr 2001 120051) . Solid lines show our full evolutionary 
calculations for a set of initial conditions representative of those 
expected in the pre-AM CVn population. Arrows on several of the 
tracks indicate the sense of evolution. The dotted lines show the 
contribution of the GW term (i.e. M = 0) in equation for the 
same M\ M2 i as the full calculations. The three lines have, from 
bottom to top'in Mq: (Mi,i, Af 2 ,i) = (0.4,0.1), (0.625,0.2), and 
(1.025,0.3). The different solid lines starting on each dotted line 
differ in initial entropy (and hence i?2,i)- The diamonds with error 
bars (only visible in J1914) show the observed locations of v for 
J0806 and J1914. 

Hz s- 1 :lStrohmayerll200l and J1914 (7 x 10~ 18 Hz s _1 ; 
lStrohmaverll2004[) . The set of initial conditions expected 
for double-degenerate channel AM CVn progenitors nat- 
urally produce secular evolution consistent with the obser- 
vations of both systems. While J0806 is consistent with 
a system whose v is dominated by GW losses, J1914 is 
only so for very low values of M\ t i, Mzy. at M = 0, if 
Mi.i > 0.15, M 2 .i < 0.1 is required. The combination 
of these lo w masses are not plau sible outcomes of binary 
evolution ijNelemans et al.ll2001|) . Thus, regardless of X- 
ray production model, a significant negative contribution 
to v (e.g., due to mass transfer or properly oriented spin- 
orbit coupling, as in the UI model) is required for J1914. 
In an accretion model, this means the M in J1914 must 
be large. 

While for Rl > R2, M g rows exponentially with 
Ai? = R L - R 2 (|Ritterl 119881). this growth rate slows 
considerably once R L < R<> (|Kolb fc Ritterl[l990j) . The 
donors also initially have non-degenerate outer layers 
that, in most cases, are predominantly radiative, pro- 
ducing contraction upon mass loss during the M turn-on 
phase. Sub-photosphcric mass transfer and donor con- 
traction both increase the duration of the turn-on phase 
beyond the previous estimate of Tj\ = M/M ^1 — 10 yrs 
made for J0806 at M = (jWillems fc Kalogeral200'5l) . In 
our full calculations, systems with P or b < 10 min reach 
i> = after 10 3 — 10 6 yrs. Thus, AM CVn binaries re- 
main in the M turn-on phase significantly longer than 
previously estimated. 

3. ACCRETION MODEL IMPLICATIONS 




Fig. 2. — The predicted v versus L a cc from our full evolutionary 
calculations for systems with the measured v and P or b of RX J0806. 
Symbol shapes indicate the initial donor mass, Mi,i as indicated 
in the Figure's key. The solid lines are included to help discern 
lines of constant M2 i. Filled symbols indicated systems that have 
always accreted via direct-impact. 



We have shown that a secular v > docs not pre- 
clude accretion driven models for J0806 and J1914. We 
now discuss the implications of accretion models for these 
systems. We first determine M^;, M 2 - l: and R 2 \ leading 
to systems consistent with the observed P or b and v in 
J08086 and J1914 and use the successful initial condi- 
tions to determine expected ranges for current system 
properties. 

We first consider v in J0806. For GW emission alone, 
this is given by 

«, , / v \ 19 / 3 ( M \ 10/ 

(2) 

where M 5/3 = (MiM 2 )/(Mi + M 2 ) 1/3 is the system's 
chirp mass. For M > 0, v < Pgr generically. For models 
passing through P or b = 321.5 and v = 3.63 x 10~ 16 Hz 
s _1 , we calculate v and L acc = M(0li —4>Ri)j where </>li; 

are the potential at the inner Lagrange point and 
surface of the accretor ijHan fc Webbinklll999lL We plot 
these predictions in Figure The lines connect mod- 
els with equal M 2j i and larger L acc values correspond to 
larger i? 2 .i- For low M, larger R 2l requires larger M to 
satisfy the v constraint. As i?2,i increases, M contribu- 
tions to v become important, reducing i) and eventually 
driving it negative. 

For J0806, the observed X-ray luminosi ty is Lx ~ 
1 — 5 x 10 32 erg s _1 (for a distance of 500 pc: llsrael et alJ 
ll999ll2003Tl . Given uncertainties in the source distance 
and the conversion efficienct of accretion energy into soft 
X-rays, we take as a rough upper limit for L acc 10 34 
erg s _1 . This places J0806 in the regime where M contri- 
butions to v arc of order or less than GW contributions. 
We predict v = 1 — 1.5 X 10~ 28 Hz s~ 2 . Given a phase 
measurement accuracy of ~ 0.01, a total timing measure- 
ment baseline of 20 yrs should be sufficient to constrain v 



Fig. 3. — The relation between L aC c at the point each model 
reaches the P or b of J0806 (lower panel) or J1914 (upper panel) and 
that model's initial P or b for models with i> equal to the observed 
value is each system. The symbols have the same meaning as in 
Figure|2]with the addition that triangles show models with M2A = 
0.1-Mq. The much narrower P or b.i range in J1914 at fixed Mi,\ is 
due to the steep dependence of ii on P or b once M is large enough 
to compete with GW emission in setting u (see Figure HI 



at this level (Strohmaycr 2006, private communication). 
Thus, this prediction should be testable within 5-10 yrs. 

We indicate by solid symbols in Figure |5] systems in 
which accretion is via direct impact. Systems which 
form a disk (open symbols) are unlikely models for J0806 
since it is then unclear how the 100% modulation of 
the X-ray light would be obtained. The direct-impact 
constraint requires M 2 ,i > 0.20M Q in J0806. We also 
check if the advection of or bital angular mom entum onto 
the accretor is significant ijMarsh et alJ l20041. We com- 
pare the GW j/J to + q)r h M 2 /M 2l where r h is 
the effective radius at which an orbit has the same spe- 
cific angular momentum as the transferred m atter (see 
iVerbunt fc RappaportllT988tlMarsh et al.ll2004^ and find 
that for L acc < 10 34 erg s _1 advection contributics < 1% 
to the J evolution. Finally for the direct-impact models, 
« 10 2 — 10 4 yrs (increasing with L acc ). 

For J 19 14, the small v requires a large M > 3 x 
1O _9 M0 yr _1 . At its longer P or b, the phase space leading 
to direct-impact systems is small: we find that M 2i \ w 
0.15 — O.2OM0 is required and that f?2,i must be close 
to the smallest value producing systems consistent with 
J1914. This M 2 _\ range is that most li kely produced by 
binary evolution ijNelemans et al.ll2001(l , so a fine-tuning 
problem is somewhat alleviated. The potential contri- 
bution of J advection onto the accretor is significant in 
J1914: (1/(1 + q)r h M 2 /M 2 J/J) « 0.4 - 0.6. Given 
the large uncertainty in tidal synchronizati on rates that 
would feed this J back into the orbit (see iMarsh et alJ 
12004 and references therein), we leave this point for fu- 
ture investigation. Our direct-impact accreting models 
have L< acc 10 34 - 10 35 and r A y « 5 x 10 5 yrs, consis- 
tent with the estimates fo r L x ~ 10 33 - 10 35 erg s~ 1 
ijRamsav et alJl2005l 120061) . 
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With the assumption that any J advection is negligi- 
ble or is cancelled by sufficiently strong tidal coupling, 
we can place constraints on initial conditions producing 
systems consistent with J0806 and J1914. In Figure [3J 
we show our constraints on P rb,i vs. L acc . 

4. SUMMARY & DISCUSSION 

We have shown that the measured v in J0806 and 
J1914 are consistent with the secular v expected for ac- 
cretion driven systems leading to AM CVn binaries. As 
such, the v measurements do not provide a means of 
distinguishing between accretion and UI models posited 
for these systems. We have shown that accretion models 
consistent with observations give r fl/ 10 2 — 10 4 yrs in 
J0806 and t,> ~ 5 x 10 5 yrs in J1914. For comparison, 
iDall'Osso et al.l l)2006f) estimate that the UI timescale for 
J1914 is « 10 4 - 3 x 10 5 yrs. We also predicted that ac- 
cretion models give v « 1 — 1.5 X 10~ 28 Hz s~ 2 for J0806; 
another 5-10 yrs baseline of timing J0806 will be suffi- 
cient to test this prediction. Quantitative predictions of 
v for the UI model are needed to determine if v will be 
a useful diagnostic tool. Our calculations also constrain 
each system's initial conditions and we relate the current 
iacc to P nr h ; in both J0806 and J 1914. 

Rccenltv. ID'Antona et"aD (|2006l) proposed that J0806 
could be modeled with a degenerate He WD donor with 
a thick H atmosphere supporting p — p burning. In their 
model, P or b decreases until the H layer is removed, expos- 
ing the degenerate He WD. Our pure He donor models 
shows that a thick H layer is not necessary to explain a 
phase where v > during early-contact AM CVn evolu- 
tion. The combination of sub-photospheric mass transfer 
and spread in initial donor degeneracy are sufficient to 
produce either J0806 and J1914. Indeed, J1914 is natu- 
rally accomodated within t he AM CVn initial cond itions 
of our model, but not bv ID'Antona et all l)2f)0(f) . prc- 
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sumably becuase they assume a fully degenerate donor 
interior. We agree, however, that the presence of H in 
J0806 or J1914 is not unexpected given the uncertainty 
in the amount of H left on the donor after the CE phase. 
Pure He models are a limiting case that could be altered 
to include an arbitrarily thick H layer, changing our re- 
sults quantitatively. As such, a detection of H in J0806 
or J1914 does not rule out accretion m odels nor on its 
own d iscriminate between our and the ID'Antona et alJ 
l)2006|) models. 

Future work on these systems must seek to diagnose 
the source of their X-ray light. An avenue for this would 
be consideration of the optical output of the systems 
where we would expect significant differences between 
the UI and accretion models. Assessment of potential 
fine-tuning problems for accretion models in J1914 is 
also needed. This would be best performed in a general 
study addressing the detection probability for early con- 
tact AM C Vns. Semi-degenera te systems are formed less 
frequently l)Delove et al.l 120051) . but evolve more slowly 
at longer P or b, so the a prioi expectations are unclear. 
This would also constrain theoretical expectations for 
the intial conditions leading to AM CVn binaries. The 
usefulness of such studies, however, hinges on the de- 
termination that J0806 and J1914 are indeed direct- 
impact accretors, not unipolar inductors. The commu- 
nity should continue efforts toward distinguishing be- 
tween these models for these enigmatic systems. 
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